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|7j " Abstract 

Context. Observations of SNRs in X-ray and y-ray bands promise to contribute with important information in our understanding of 
the kinematics of charged particles and magnetic fields in the vicinity of strong non-relativistic shocks and, therefore, on the nature 
of galactic cosmic rays. The accurate analysis of SNRs images collected in different energy bands requires the support of theoretical 
modeling of synchrotron and inverse Compton emission from SNRs. 

Aims. We develop a numerical code (remlight) to synthesize, from MHD simulations, the synchrotron radio, X-ray and inverse 
Compton y-ray emission originating from SNRs expanding in non-uniform interstellar medium (ISM) and/or non-uniform interstellar 
magnetic field (ISMF). As a first application, the code is used to investigate the effects of non-uniform ISMF on the SNR morphology 
in the non- thermal X-ray and y-ray bands. 
Methods. We perform 3D MHD simulations of a spherical SNR shock expanding through a magnetized ISM with a gradient of 
. ambient magnetic field strength. The model includes an approximate treatment of upstream magnetic field amplification and the effect 

of shock modification due to back reaction of accelerated cosmic rays, assuming both effects to be isotropic. From the simulations, 
we synthesize the synchrotron radio, X-ray and inverse Compton y-ray emission with the synthesis code remlight, making different 
assumptions about the details of acceleration and injection of relativistic electrons. 

Results. A gradient of the ambient magnetic field strength induces asymmetric morphologies in radio, hard X-ray and y-ray bands 
independently from the model of electron injection if the gradient has a component perpendicular to the line-of-sight (LoS). The 
degree of asymmetry of the remnant morphology depends on the details of the electron injection and acceleration and is different in 
the radio, hard X-ray, and y-ray bands. In general, the non-thermal X-ray morphology is the most sensitive to the gradient, showing 
the highest degree of asymmetry. The IC y-ray emission is weakly sensitive to the non-uniform ISMF, the degree of asymmetry of the 
remnant morphology being the lowest in this band. 

Key words. Magnetohydrodynamics (MHD) - Radiation mechanisms: non-thermal - Shock waves - ISM: supernova remnants - 
Gamma rays: ISM - X-rays: ISM 

1- Introduction | 2001f lAlbert et al] l2007b RX J17 13.7-3946. iMuraishi et alJ 

l2000t lEnomoto et al.l l2002t lAharonian et al.l l2007aT 



•£] ■ It is largely accepted m the literature that the hard X-ray emis- Berezhko & V 6lk 2008; RX J08 52.0-4622. iKatagiri et al.l 

K>S sion detected m many young shell-type supernova remnants [^OoIlAharonian et al. 2005, 2007b; Enomoto et al. 20o3lnd 

53 . (SNRs) is synchrotron emission from electrons accelerated to RCW J^ 0VVQ & Lemoine-Goumard 2008) allows both for 

energies of tens of TeV (e.g. | Koyama et al. || 1995 | ) by means of leptonic and hadronic origin of VHE y-rays. 
the diffusive shock acceleration process. In addition, inverse 

Compton (IC) collisions of these high energy electrons with Recently, the study of SNRs as particle accelerators has re- 
low energy photons from the ambient radiation field (e.g. the ceived a strong impulse thanks to the new y-ray observations of 
cosmic microwave background; hereafter CMB) are expected, SNRs with the instruments of the High Energy Spectroscopic 
thus leading to very-high-energy (VHE; > 100 GeV) y-ray System (HESS), the Major Atmospheric Gamma-ray Imaging 
emission too. In regions of high mass density, ions in the shell Cherenkov (MAGIC) experiments and the Fermi Gamma-ray 
are likely to have been accelerated to similar energies, and Space Telescope. The analysis of multi-wavelength observa- 
y-rays may be due to neutral pion decay from proton-proton tions, from radio, to hard X-rays, to y-rays, promises to increase 
interactions. The nature of the TeV emission, therefore, stands our understanding of the kinematics of charged particles and 
on the combination of X-ray synchrotron emitting electrons magnetic fields in the vicinity of strong non-relativistic shocks 
and very energetic ions; it is not clear, at the present time, and of the possible role of SNRs for the origin of galactic cosmic 
which one between the two mainly contributes to the de- rays (CRs). In this context, a very important source of informa- 
tected VHE y-rays. The spectral analysis of multi-wavelength tion could be the distribution of surface brightness observed in 
data of several shell-type SNRs (Cas A, lAharonian et al.l SNRs in several bands. For instance, the properties of the bright- 

ness distribution have been crucial in our understanding of the 

Send offprint requests to: S. Orlando, acceleration and injection of relativistic electrons thr ough SNR 

e-mail: orlando@astropa.inaf.it shocks (e.g. the criterion of iRothenflug et al] 12004] versus the 
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azimuthal profile comparison in iPetruk et al.l l2009c). Recently, 
Petruk et al. (2010a, submitted to MNRAS) have compared syn- 
thetic distributions of surface brightness predicted for SN 1006 
in different bands with observations, deriving important observa- 
tional constraints on the modeling of SN 1006. Also, the corre- 
lations of brig htness in radio, X-r ays and y-rays claimed in RX 
J1713.7-3946 (lAcero et alJ l2009b and some others SNRs (e.g. 
Aharonian et al . 2006) could be considered to favor electrons as 
being responsible for VHE y-rays in these SNRs. 

The analysis of the brightness distributions observed in 
different energy bands needs to be supported by an accurate 
comparison of the observed distributions with those predicted 
by detailed MHD models. For insta nce, a number of SNRs 
present a bilateral structure (BSNRs; | Kesteven & Caswellll 19871: 
iFulbright & Revnoldslfl990l: lGaenslerlll998h and it is not clear 
how to "translate" this 2D information into a 3D morphology 
of the emission over the SNR shell. The two competing cases 
traditionally invoked are either equatorial-belt or polar-caps and 
they are related to the model of injection of relativistic electrons 
(isotropic or quasi-perpendicular in the former case and quasi- 
parallel in the latter case). Establishing the 3D morphology of 
SNRs, therefore, may give some important hints on the acceler- 
ation theory. 

The distribution of surface brightness of synchrotron emis- 
sion in SNRs expanding through a uniform interstellar medium 
(ISM) and uniform interstellar magnetic field (ISMF) has been 
extensivel y investigated, through nume rical modeling , in both 
the radio (iFulbright & Reynoldsl Il990h and X-ray (iReynoldsl 
1 19981 12004]) bands; in particular, the dependence of the bright- 
ness distributions on the efficiency of the acceleration process 
has been explored, considering different injection models. First 
IC y-ray ma ps of SNRs in a unifo rm ISM and ISMF have been 
presented by IPetruk et aD (l2009al) who investigated the proper- 
ties of brightness distributions in VHE y-rays as compared to 
the distributions in the radio band (see also Petruk et al. 2010b, 
submitted to MNRAS). Recently. IPetruk etal.1 (l2009bl) proposed 
a method to predict IC y-ray images of SNRs starting from 
observed synchrotro n radio maps and spatially resolved X-ray 
spectral analysis (e.g. lMiceli et al.ll2009l) of SNRs. 

In a previous work (lOrlando et al.ll2007i hereafter Paper I), 
we have investigated the origin of asymmetries in the radio mor- 
phology of B SNRs through a model of a SNR expanding through 
either a non-uniform ISM or a non-uniform ISMF. In this paper, 
we extend our analysis to the non-thermal X-ray and IC y-ray 
emission. In particular, we develop a numerical code (remlight) 
to synthesize the synchrotron radio, X-ray, and IC y-ray emis- 
sion from 3D MHD simulations; then we couple the synthesis 
code with the MHD model introduced in Paper I (extended to 
include a simple treatment of upstream magnetic field amplifi- 
cation and the effect of shock modification due to back reaction 
of accelerated CRs) and investigate the effects of a non-uniform 
ISMF on the morphology of the remnant in the hard X-ray and y- 
ray bands. Though remnants of Type la supernovae are expected 
to expand in an almost uniform ISMF, here we show that even a 
very small gradient of the ISMF can influence significantly the 
non-thermal remnant morphology. 

In Sect. [2] we describe the MHD model and the numerical 
setup; in Sect. [3] we describe the computation of synchrotron X- 
ray and IC y-ray emission, including model of relativistic elec- 
tron behavior; in Sect. |4] we discuss the results; and, finally, in 
Sect. [5] we draw our conclusions. 



2. MHD modeling and numerical setup 

We adopt the MHD model introduced in Paper I, describing 
the propagation of a SNR shock through a magnetized ambi- 
ent medium. The shock propagation is modeled by numerically 
solving the time-dependent ideal MHD equations of mass, mo- 
mentum, and energy conservation in a 3D Cartesian coordinate 
system (x, y, z) (see Paper I for details). The model does not in- 
clude consistently the effects on shock dynamics due to back- 
reaction of accelerated CRs. However, we approach the effect 
of shock modification by considering different values of the adi- 
abatic index y which is expected to drop from the value of an 
ideal monoatomic gas; in particular, we consider here the cases 
of y = 5/3 (for an ideal monoatomic gas), y = 4/3 (for a gas 
dominated by relativistic particles), and y = 1.1 (for large en- 
ergy drain from the shock region due to the escape of high en- 
ergy CRs). In addition, we account for upstream magnetic field 
amplification due to back reaction of accelerated protons, by am- 
plifying the (pre- shock) ambient magnetic field in the neighbor- 
hoods of the remnant. These effects, namely shock modification 
and magnetic field amplification, might depend on the obliq- 
uity angle (i.e. the angle between the external magnetic field 
and the normal to the shock); for instance, they could follow 
the same dependence of the injection efficiency (i.e. the frac- 
tion of accelerated electrons). However, at the present time, the 
dependence of these effects on the obliquity angle is not un- 
derstood and, therefore, we consider here the simplest case of 
isotropic shock modification and magnetic field amplification 
(i.e. no additional obliquity-dependent magnetic-field amplifica- 
tion has been assumed). On the other hand, such a choice makes 
easier our analysis of non-thermal images which are already in- 
fluenced by the obliquity dependence of other processes (e.g. in- 
jection efficiency, magnetic field compression, maximum energy 
of electrons; s ee Sect. [3]). The sim ulations are performed using 
the flash code (iFryxell et al.ll2000h . an adaptive mesh refinement 
multiphysics code for astrophysical plasmas. 

As initial conditions, we adopt parameters appropriate to re- 
produce the SNR SN 1006 after 1000 yr of evolution: we assume 
an initial spherical remnant with radius ro sn r = 0.5 pc, originat- 
ing from a progenitor star with mass of 1.4 M sun , and propa- 
gating through an unperturbed magneto- static medium. The ini- 
tial total energy Eq is set to a value leading to a remnant ra- 
dius r snr « 9 pc at t = 1000 yr (Eo ranges between « 1.3 and 
1.8 x 10 51 erg, depending on y) and is partitioned so that most of 
the SN energy is kinetic energy. The remnant expands through 
an homogeneous isothermal medium with particle number den- 
sity n = 0.05 cm -3 and temperature T = 10 4 K. We consider 
three different configurations of the unperturbed ambient mag- 
netic field: 1) a uniform ambient magnetic field (runs Unif-gl, 
Unif-g2, and Unif-g3); 2) a gradient of ambient magnetic field 
strength perpendicular to the average magnetic field (runs Grad- 
BZ-gl, Grad-BZ-g2, and Grad-BZ-g3); and 3) a gradient of am- 
bient magnetic field strength aligned with the average magnetic 
field (runs Grad-BX-gl, Grad-BX-g2, Grad-BX-g3). 

In the case of a uniform ISMF, we assume that the field is ori- 
ented parallel to the x axis. In the other two cases, the ambient 
magnetic field is assumed to be dipolaiQ. The dipole is oriented 
parallel to the x axis and located either on the z axis (x = y = 0) 
at z = -100 pc (Grad-BZ-gl, Grad-BZ-g2, Grad-BZ-g3) or on 
the x axis (y = z = 0) at x = -100 pc (Grad-BX-gl, Grad-BX- 
g2, Grad-BX-g3). In all the cases, we assume the magnetic field 
strength to be Bq = 30 yuG at the center of the SN explosion 

1 This idealized situation is adopted here mainly to ensure magneto- 
staticity of the non-uniform field. 
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Table 1. Relevant initial parameters of the simulations. 







r 


[10 51 erg] 


magnetic field 
configuration 


(x,y, Z y 
pc 


Unif-gl 




5/3 


1.30 


uniform 




Unif-g2 




4/3 


1.54 


uniform 




Unif-g3 




1.1 


1.81 


uniform 




Grad-BZ-^ 




5/3 


1.30 


z-strat. 


(0, 0,-100) 


Grad-BZ-| 




4/3 


1.54 


z-strat. 


(0,0,-100) 


Grad-BZ-| 




1.1 


1.81 


z-strat. 


(0, 0,-100) 


Grad-BX- 


gl 


5/3 


1.30 


x-strat. 


(-100, 0, 0) 


Grad-BX- 


§2 


4/3 


1.54 


x-strat. 


(-100, 0, 0) 


Grad-BX- 


§3 


1.1 


1.81 


x-strat. 


(-100, 0, 0) 



a Initial energy of the explosion. b Coordinates of the magnetic dipole 
moment. 



(x = y = z = 0), roughly an order of magnitude higher than the 
galactic ISMF expected at the location of SN 1006. This high 
value has been chosen to mimic the effects of upstream mag- 
netic field amplification (see discussion above). In such a way, 
the post-shock ambient magnetic field is expected to be of the 
order of a few hundred yuG as deduced from observations. In the 
configurations with non-uniform ISMF, the field strength varies 
by a factor ~ 6 over 60 pc: either in the direction perpendic- 
ular to the average ambient field (B) (Grad-BZ-gl, Grad-BZ- 
g2, Grad-BZ-g3), or parallel to <B> (Grad-BX-gl, Grad-BX-g2, 
Grad-BX-g3). We follow the expansion of the remnant for 1000 
yr. Table [T] summarizes the physical parameters characterizing 
the simulations considered here. 

The SN explosion is at the center (x, y, z) = (0, 0, 0) of the 
computational domain which extends between -10 and 10 pc in 
all directions. At the coarsest resolution, t he adaptive mesh algo- 
rithm used in the flash code (paramesh; iMacNeice et al.ll2QQQb 
uniformly covers the 3D computational domain with a mesh of 
8 3 blocks, each with 8 3 cells. We allow for 5 additional nested 
levels of refinement during the first 100 yr of evolution with reso- 
lution increasing twice at each refinement level; then the number 
of nested levels progressively decreases down to 2 at t = 1000 
yr as the remnant radius increases following the expansion of the 
remnant throug h the mag netized medium. The refinement crite- 
rion adopted (Lohner 1987) follows the changes in density and 
temperature. This grid configuration yields an effective resolu- 
tion of « 0.0098 pc at the finest level during the first 100 yr 
of evolution (when the radius of the remnant was < 2 pc) and 
« 0.078 pc at the end of the simulation, corresponding to an 
equivalent uniform mesh of 204 8 3 and 25 6 3 grid points, respec- 
tively. We assume zero-gradient conditions at all boundaries. 

The model does not include the radiative cooling, describing 
only the free and adiabatic expansion phases of the remnant. The 
tran sition time from adiabatic to radi ative phase for a SNR is 
(e.g. lBlondinet al. 1998; Petruk 2005) 



t ti = 2.84 x 10 z 



c4/17 -9/17 

E' n. ' yr , 

51 ism J ' 



(1) 



where E51 = Eo/(l0 51 erg) and n{ sm is the particle number den- 
sity of the ISM. In our set of simulations, t it > 10 5 yr which is 
much larger than the time covered by our simulations. Our mod- 
eled SNRs therefore never reach the radiative phase. On the other 
hand, here we aim at describing young SNRs (i.e. those that are 
still in the adiabatic expansion phase) from which non-thermal 
X-ray and y-ray emission is commonly detected. 



3. Synchrotron X-ray and inverse-Compton y-ray 
emission (REMLIGHT) 

From the model results, we synthesize synchrotron rad io, X-ray, 
and IC y-ray emission, by generalizing the approach of lReyn olds 
(1998) to cases of non-uniform ISM and/or non-uniform ISMF. 
In Paper I, we have already discussed the synthesis of syn- 
chrotron radio emission and we refer the reader to that paper for 
the details of calculation. Here we discuss the synthesis of X-ray 
and IC y-ray emission as it is implemented in the synthesis code 

REMLIGHT. 

One could assume that the synchrotron X-ray or IC y-ray ra- 
diation is due to relativistic electrons distributed with an energy 
spec trum N(E) = KE~ S exp(-E/E mSiX ) electrons cm -3 erg s" 1 
(e.g. iGaisser et al.ll 199 8), where E is the electron energy, N(E) 
is the number of electrons per unit volume with arbitrary direc- 
tions of motion and with energies in the interval [E, E + dE], K 
is the normalization of the electron distribution, s the power law 
index, and £ max the maximum energy of electrons accelerated by 
the shoclQ. Nevertheless, some observations suggest that the cut- 
ofT co uld be broader than pure exponent (e.g.lEllison et al.l l2000. 
l200ll:IUchivama et al.l2003l:lLazendic et alJ2004h . Therefore we 
assume that the energy spectrum of electrons is given by 



N(E) = KE~ S exp 



(2) 



where a < 1 is the parameter regulating the broadening of the 
high-energy end of electron spectrurr0. 

The volume emissivity due to synchrotron or IC radiation 
can be expressed as 



f 

Jo 



N(E)A(E,e, B) dE, 



(3) 



where A(E, s, B) is the radiation power of a single electron with 
energy E, and s is the photon energy. The emissivity, i(e), de- 
pends on the magnetic field strength, B, only in the synchrotron 
emission process. We compute the surface brightness of the SNR 
at a given energy s, by integrating the emissivity i(e) at each 
point along each LoS in a raster scan (assuming that the source 
is optically thin). 

In the case of synchrotron emissivity in the X-ray band, 
z'x(£), the spectral distribution of radiation power of a single elec- 
tron with energy E in the magnetic field B is 



A x (E,e) = — F\ — 



(4) 



where s c = hv c = hci^BE 2 , h is the Plank constant, v c is the 
critical frequency, the angle between the magnetic field and 
the LoS, ji^ is either - sin0 for the case of ordered mag- 
netic field or ^ = <sin0) = n/4 for disordered magnetic field, 
c\ - 3e/(4nmlc 5 ), e and m e are the charge and mass of electron, 
respectively, c is the speed of light. The special function F(w) 



2 Note that the distribution of electrons N(E) in the case of syn- 
chrotron radio emission is expressed as N(E) = KE~ S electrons cm -3 
erg s _1 (see, for instance, Paper I). 

3 Note however that recent work by Kang& RvT] ( l20l0h suggests that 
or > 1. 
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can b e approximated as (e.g. lRybicki & Ligh tman 1985; Wallis 



F(w) = 



f2.15w 1/3 w<0.01, 
V7rw a29 exp(-w) 0.01 < w < 5 , 
^/n/2w l/2 exp(-w) w > 5 . 



(5) 



We found the above approximation quite accurate with discrep- 
ancies <4% from the exact value. In addition, JT Fds = 1.59 

while the exact value is 87r/9 V3 = 1.61. 

In the case of y-ray emissivity due to IC process, z'ic, the 
spectral distribution A\c(E,s) of radiation power of a single 
electron in a blac k-body photon field in Eq. [3] is (see also 
IPetruk et al . 2009a) 



2e 4 e c 



nfi 3 c 2 



r licirjcTJo) , 



(6) 



where Y is the Lorenz factor of electron, e c - &7cmbr, ?cmbr is 
the temperature of the cosmic microwave background radiation 
(CMBR) assumed to be Tcmbr = 2.75 K, 



6c S 



(m e c 2 ) 



,2\2 ' 



?7o 



4rm e c 2 (rm e c 2 - s) 



(7) 



and the s pecial functio n Ji C (?7c>?7o) can be accurately approxi- 
mated as (Petruk 200|) 



Iic(jic,m) ~ "^ciexp 



+2//o exp 



1/2 



> exp 



2/^o 
3Uc 



(8) 



This approximation represents Ji c (?7c, ?7o) in any regime, from 
Thomson to extreme Klein-Nishina. The approximation is ex- 
act in the Thomson limit. It restores detailed calculations with 
maximum error of 30% in the range of parameters which gives 
non-negligible contribution to emission. 

3.1. Maximum energy of electrons 



We follow the approach of lReynoldsl (Il998l) for the description 
of time evolution and surface variation of E max , generalizing 
his ap proach to cases o f non-uniform ISM and/or non-uniform 
ISMF. iReynoldsl (Il998l) considered three alternatives for time 
and spatial dependence of E max . Namely, the maximum accel- 
erated energy maybe determined: 1) by the electron radiative 
losses (due to synchrotron and IC processes), 2) by the limited 
time of acceleration (due to the finite age of the remnant) and 3) 
by properties of micro-physics when the scattering of electrons 
with E > £ max becomes less efficient and the electrons freely 
escape from the region of acceleration^. The maximum energy 
is given by 

iW <* /e*(©o) V* B A J , (9) 

4 For the escape case, it is commonly assumed that MHD waves re- 
sponsible for the scattering are much weaker above some wavelength, 
A max , and E max is approxim ately the energy of particles with that gyro- 
radius (e.g. Reynolds 1998). 



where is the obliquity angle, /e,^(© ) is a function describ- 
ing smooth variations of £ max versus obliquity, V S h is the shock 
velocity, B is the pre-shock ISMF strength, and £ = 1, 2, 3 cor- 
responds respectively to loss-limited, time-limited and escape- 
limited models of E max . The values of q and A are: q\ = 1, 
q 2 = q 3 = 0, and A\ = -1/2, A 2 = A 3 = 1 dReynolddll998h . 
Note that we assume q 2 - because 1) £ max rises quite slowly 
with time when it is determined by the finite time of accelera- 
tion (Re ynoldsll 19981) . even in the nonuniform ISM, and 2) most 
of emission rises from downstream regions close to the shock; 
therefore, the variation of £ max due to velocity variation is neg- 
ligible in the (thin) emitting region. From Eq. we express the 
surface variation of £ max as 



iW = ^max^ll /e*(®o) ^ > 



(10) 



where £max,£|| is a free parameter, representing the maximum en- 
ergy in a point p on the SNR surface where the ISMF is parallel 
to the shock normal, Tsh = V sh /V s h,ih = B /B 0j \\, V sh ,n and 
2?o,|| are the shock velocity and pre-shock ISMF strength in the 
point p, respectively. 

A detailed theoretical framework p roviding the obl iquity 
dependence of £ max was presented by Reynolds (1998) and 
is based on the prescription for diffusion by iJokipiil (Il987l) . 
However, this theory is limited to the test-particle regime, as- 
suming no magnetic field amplification. On the other hand, at 
the present time, a more general theory describing the obliquity 
dependence of £" max is still lacking. For the sake of generality, 
we adopt here some arbitrary smooth variations of £ max versus 
obliquity with the goal to see how different trends in the obliq- 
uity dependence of £ max influence the visible morphology of 
SNRs. In fact, the remnant morphology - once we are not in- 
terested in the exact comparison with observations - is mainly 
affected by the contrast C ma x = £max,||/£max,± and not by the 
exact shape of the dependence on obliquity, once the latter is 
assumed to be smooth. 

Our strategy is to consider smooth variations of £ max ver- 
sus obliquity that correspond to the loss-limited, time-limited 
and esc ape-limited mod els of E max in the theoretical frame- 
work of Reynolds (1998). Since, in general, the mechanism lim- 
iting electron acceleration (i.e. loss-limited, time-limited and 
escape-limited) may be different at different shock obliquity 
angles, we calculate the value of £ max appropriate for each 
limitation mechanism at each point, by considering E max = 
min[£ max j, £ m ax,2, £max,3] (where indexes 1, 2, 3 correspond re- 
spectively to loss-limited, time-limited and escape- limited m od- 
els). This way to compute £ max is adopted in Sect. 14.21 14.31 and 
14.6 1 w here we assume also to be in the Bohm limit (i.e. gyrofac- 
to£| 7/ = 1) in the test-particle regime. In particular, in Sect. 14.21 
we introduce a reference case for which the adopted set of pa- 
rameters (see Sect. 14.11) leads to C ma x > 1 • Note that the adopted 
parameters make this case suitable for comparison with young 
non-thermal SNRs as, for instance, SN 1 006. In addition, for the 
sake of generality, in Sect. 14.41 and 14.51 we explore the effects 
on the remnant morphology of various obliquity dependencies 
of £max, by considering also cases for which the contrast C ma x is 
< 1. 



5 The "gyrofactor" is defined as the ratio between th e mean free path , 
A\\, along the magnetic field and the gyroradius, r g (see Reynolds 1998). 
In general it is expected that the mean free path can be no less than r g , 
so that r] > 1; the equality corresponds to the Bohm limit, i.e. a level 
of turbulence leading to wave amplitudes comparable to the stationary 
magnetic field strength. 
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3.2. Post-shock evolution of the electron distribution 

As in iRevnoldsl (119981) . we assume that relativistic electrons are 
confined in the fluid elements which advect them from the region 
of acceleration. Fluid element with Lagrangian coordinate a = 
R(t{) was shocked at time t u where R is the radius of the shock. 
At that time, the electron distribution on the shock was 



N(E { , h) = K s (a, h)Er s exp 



(ii) 



where E x is the electron energy at time t{, K s is the normalization 
of the electron distribution immediately after the shock (in the 
following, index "s" refers to the immediately post-shock val- 
ues), and s is the power law index. At variance with Paper I, 
we are interested here in synchrotron X-ray and IC y-ray emis- 
sion. In this case, the evolution of the electron distribution has to 
account for energy losses of electrons due to both adiabatic ex- 
pansion and radiative losses caused by synchrotron and IC pro- 
cesses. At time t- u the energy of the electron confined in the fluid 
element with Lagrangian coordinates a = R(t{) was (cf. Eq. 26 
in lRevnoldsl[l99g 



where cr = (y + l)/(y - 1) is the shock compression ratio and r 
is the Eulerian coordinate, together with the continuity equation 
p (a)a 2 da = p(a, t)r 2 dr and the derivative 



~dE 



1 



£ad£; 



(18) 



rad 



implies that downstream 

N(E,a,t) = K{a,t)E- s 8 s ~l 
x exp 



E m ax(t, ©o) T'ia, R) £ a d£rad 



(19) 



with K(a,t) = K s (ti)6 s a f, E max (t,® Q ) is given by Eq. [H and 
taking into account Eq.[l5l Assuming that K s oc p s V S h(t)~ b , i.e. it 
varies with the shock velocity V s h(0 and, in case of non-uniform 
ISM, with the immediately post-shock value of mass density, p s , 
the downstream variation of K(a, t) is described by the relation 
(see Paper I) 



Ei 



£ad^i 



(12) 



'ad^rad 



where E is the electron energy at the present time t, £ a d is a term 
accounting for the energy losses of electrons due to adiabatic 
expansion 

\Ps(*i) / \ Ps(0 / \Po(a)J 

p is the mass density (in the following, index "o" refers to the 
pre- shock values), £ ra d is a term accounting for the radiative 
losses of electrons 



£ ra d(£, a,t) = 1 - I(a, t) 



(14) 



637/(^ ff s || t) erg is the fiducial energy at parallel shock 



(Reynolds 1998), 

^eff s || — ^CMB 

is an "effective" magnetic 
field at parallel shock accounting energy losses due to IC scat- 
terings on the photons of CMB, Bqmb = 3.27 jiG is the magnetic 
field strength with energy density equal to that in the CMB, t is 
the time, and I(a, t) is an integral independent of E which is cal- 
culated with the approach described in Appendix|Al The electron 
energy losses in a given fluid element are mainly due to radiative 
losses if Ef < E max and to adiabatic expansion if Ef > E max . 

At time t{, the shock was able to accelerate an electron 
confined in the fluid element with Lagrangian coordinate a to 
Emaxik)- From Eq.[9j we derive that 

£maxfe) = / E (0ofe)) ( Vsm \ q (B (a)\ A 

£max(0 /e(®o«) UshW/ \B (R)) ^ ' 

The ratio, V s h( fi)/Vsh(0» m ^Y be expre ssed through pressure, P, 
and density, p dHnatvk & Petruk| [l999) 



VM) = ( P(a,t) \ l/2 (p (g)X 

y sh (0 \ p s (t) J \ Po (R)J 



^- [)/2 ( P (a,t))-^ 2 

Ps(0 



(16) 



The conservation law for the number of particles per unit 
volume per unit energy interval 



N(E,a,t) = N(Ei,a,tO 



a 2 da dE{ 
err 2 dr dE ' 



(17) 



K(a,t) f g (&o(ti)) ( P(a,t) 



K s (R,t) 



f g (® (t)) \Ps(RJ) 

Po(R) 



-b/2 



P(a, t) V 
Ps(R,t)J 



by/2+(s+2)/3 



(20) 



where f g (S ) is the obliquity dependence of the injection effi- 
ciency g (the fraction of accelerated electrons). Again, in the lack 
of theoretical dependence of the injection efficiency on obliq- 
uity in case of efficient acceleration, we consider a number of 
smooth functions for f g (S ), exploring either increasing or de- 
creasing injection. In fact, such an approach, which considers 
different contra sts f c \\/f c ± e i ther > 1 or < 1 is realized in our 
previous paper Petruk etal. (2009a); there, it was shown how 
the change in the injection contrast influences non-thermal im- 
ages of SNRs in uniform ISM and uniform ISMF. In the present 
paper, for simplicity and in order to see the effects from nonuni- 
formity of ISMF, we limit our considera t ions to three contrasts of 
injection. Namely, following IRevnoldsl ([1998), we consider the 
following models: quasi-parallel (J g (S ) = cos 2 @ s ), isotropic 
(f g (®o) = 1), and quasi-perpendicular (f g (® ) = sin 2 @ s ), where 
@ s is related to @ through the expression cos @ s = cr~ l cos @ - 
The first injection model leads to a three-dimensional polar-caps 
structure of the remnant, whereas the latter two produce a three- 
dimensional equatorial-belt structure of the remnant. 

Radiative losses of electrons E oc E 2 are mostly effective 
in modification of the distribution N(E, a, t) around E ~ £ max 
(Reynolds 1998). This may be noted in Eq. [T9j The variation 
of the energy distribution N(E,a,t) of electrons with energy 
E E mSLX (in this case also E «c Ef, leading to 6 ra d — > 
1), i.e. electrons with negligible radiative losses, is given by 
N(E, a, t)/N s (E, R, t) = K(a, t)/K s (R, t), where N S (E, R, t) is the 
energy distribution of electrons immediately after the shock. 
This expression does not depend on energy E and, in fact, we 
have used this expression in Paper I for investigation of prop- 
erties of surface brightness distribution of SNR emitting ra- 
dio frequencies. In contrast, the modification of the distribution 
N(E 9 a, t) due to effective electron radiation is given by the two 
last multipliers in Eq.[l9] The radiative losses of electrons there- 
fore are important for the surface brightness distribution of SNR 
in X-ray and y-rays. 
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Figure 1. 3D rendering of the mass density at t = 1000 yr for a 
remnant expanding through a uniform ISMF and for three cases 
of y = 5/3, 4/3, 1.1 (runs Unif-gl, Unif-g2, and Unif-g3, see 
Table©. 



4. Synchrotron and IC images of SNRs expanding 
through non-uniform ISMF 

The evolution of the remnant expanding through the non- 
uniform ISMF has been described in Paper I where the reader 
is referred to for more details. Figure [T] shows the 3D render- 
ing of the mass density at t = 1000 yr in the three cases of y 
considered for uniform ISMF. The main effect of y on the shock 
dynamics is to change its compression ratio and the distance of 
the contact discontinuity from the blast wave position; no depen- 
dence on the obliquity angle is present, y being uniform in each 
simulation. The value of y is expected therefore to influence the 
absolute values of emission in the radio, X-ray and y-ray bands 
but not the large scale morphology of the remnant to which this 
paper is focused on. In the following, we first discuss the ef- 
fects of non-uniform ISMF on the synchrotron and IC emission 
adopting, as reference, the case with y = 5/3, allowing the direct 
comparison of ou r results with those available in the literature; 
then in Sect. 14.31 we discuss the effect of y on the morphology 
of the non-thermal emission. 

In all the synthetic images presented below, we introduce 
the procedure of magnetic field disordering (with randomly ori- 
ented magnetic field vector in each point) downstream of the 
shock (see Paper I), according to observati ons showing a low 
degree of polarization (10-15%; e.g. Tycho, iDickel et all 1 19911 
SN 1006. iRevnolds & Gilmorel 119931) . Note that the obliquity 
angle @ s is derived from @ through the conversion formula 
cos@ s = cr~ l cos@ given in Sect. 13.21 and, therefore, does not 
take into account the m agnetic field disordering; as discussed by 
Fulbrigh t & Reynolds! (Il990h . this corresponds to the assump- 
tion that the disordering process takes place over a longer time- 
scale than the electron injection which occurs in the close prox- 
imity of the shock. 

In all the simulations, we assume the (average) unperturbed 
ISMF (B) oriented along the x axis. In the two magnetic field 





A v 






vb\ 


0VB 




LoS 




0B 











Figure 2. Relevant angles describing the orientation of the ISMF 
and of the gradient of ISMF strength with respect to the observer: 
0b is the angle between the (average) unperturbed ISMF and 
the LoS, and 0vb is the angle between the gradient of the ISMF 
strength and the vertical line passing through the center of the 
remnant Ay. 



configurations explored in this paper, the gradient of ISMF 
strength is either normal (runs Grad-BZ-gl, Grad-BZ-g2, Grad- 
BZ-g3; V|(B)| along z) or aligned (runs Grad-BX-gl, Grad-BX- 
g2, Grad-BX-g3; V|<B)| along x) to (B). Since we analyze the 
remnant morphology as it would be observed from different 
points of view, we define two angles to describe the orientation 
of (B) and V|(B)| in the space (see Fig. 0: 0b is the angle be- 
tween (B) and the LoS, and 0vb is the angle between V|(B)| 
and the normal to the ISMF in the plane of the sky (axis A v in 
Fig.O. The first angle is the aspect angle commonly used in the 
literature. The definition of the second angle allows us to explore 
the remnant morphology for various aspect angles and for fixed 
0vb, V |(B) | lying on a cone with angle 0vb (see Fig.0. In cases 
in which the gradient VB is aligned with the average ISMF (runs 
Grad-BX-gl, Grad-BX-g2, Grad-BX-g3) V b = 90° by defini- 
tion. In Grad-BZ models, the angle between (B) and V|(B)| is 
always 90°. In the following, the images are calculated for vari- 
ous values of the angles defined above and with a resolution of 
256 x 256 pixels. 

4.1. Parameter space 

The prescriptions for the electron energy distribution at any point 
inside the remnant and for the synthesis of synchrotron and 
IC emission discussed in Sect. [3] are characterized by several 
parameters regulating the energy spectrum of relativistic elec- 
trons, the injection efficiency, the time and spatial dependence 
of ^max> etc.. In the following, we limit the model parameter 
space through some assumptions that allow us to fix some of the 
parameters. 

In particular, we assume that the power law index in Eq. [2] 
is s = 2, as sugg ested by man y observations of BSNRs (e.g. 
for SN 1006: lMiceli et al.ll2009h . In the test-particle regime, the 
index should be related to the shock compression ratio through 
s = (2 + cr)/(cr - 1) with cr = (y + l)/(y - 1). In case of efficient 
shock acceleration, the electron energy distribution is curved. 
The value of s that should be used in an approximation like Eq.[2] 
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is rather a m ean radio-to-X-ra y spectral index whose value is 
around 2 (e.g. I Allen et al.ll2008l) . independent of the local slopes 
of the electron spectrum. As for the curvature of the spectrum 
around E max , we assume that a = 0.5 in Eq.[2] In fact, this could 
be the case in SN 1006 and G347.3-0.5 where a number of mod- 
els su ggest a ~ 0.5 jEUison et all 12000. 200 it lUchiyama et"aD 



els su ggest a « u.;> qbiiisor 
2003: lLazendic et alJ l2004). 



The maximum energy at parallel shock £ max ,|| is a free pa- 
rameter in Eq.[TOlthat we assume to be £ m ax,n = 26 TeV in most 
of our calculations. This parameter has to be compared with the 
fiducial energy at parallel shock defined as Ef t \\ = 637 /(B^ s ^ t) 
erg (lReynoldslll998l) ; for the cases considered here, we have 
E u = 14 TeV in models with uniform ISMF, E u = 12 TeV 
in Grad-BZ models and Ef t \\ = 2 TeV in Grad-BX modelfl In 
all these cases, therefore, Ef < E mSLX for a significant portion 
of the remnant and the electron energy loss es are m ainly due to 
radiative losses (see discussion in Sect. 13.2b . In Sect. 14.61 we in- 
vestigate the dependence of the non-thermal emission on £ m ax,ih 
by exploring cases for which E mSLX < Ef and the electron energy 
losses are mainly due to adiabatic expansion. 

The parameter b in Eq.[20|is a constant and determines how 
the injection efficiency depends on the shock properties; we as- 
sumed that K s oc p s V S h(t)~ b (see Sect. 13. 21 and Paper I). On theo- 
retical grounds b might be expected to be negative, reflecting an 
expectation that injection efficiency may behave in a way similar 
to acceleration e fficiency: s tronger shocks might inject particles 
more effectively. iReynoldsl ([1998) considered three empirical al- 
ternatives for b as a free parameter, namely, b = 0,-1, -2. In 
particular, b = -lis commonly assumed in many areas of astro- 
physics such as gamma-ray bu rsts and prompt rad i o and X-ray 
emission from SNe. However iBandiera & Petrukl (l2010h have 
shown that models preferring a constant fraction of the shock 
energy to be transferred into CRs (i.e. b = -2) are rejected by 
statistical analysis of two SNR samples. In addition, Petruk et 
al. (2010a, submitted to MNRAS) compared their model results 
with experimental data of the remnant SN 1006 and found that 
b has a value between and -I. Petruk et al. (2010b, submit- 
ted to MNRAS) showed that the smaller b, the thicker the radial 
profiles of the surface brightness in all bands; an effect mostly 
prominent in radio band. Since no effect on the pattern of asym- 
metries induced by a non-uniform ISMF is expected (see Paper 
I), we assume b - in all our calculations, this being the most 
neutral case. 

The synthe tic images are expected to depend on the remnant 
age (iReyno lds 1998). To reduce further the number of model pa- 
rameters, we focus here on remnant 1 kyr old, as in the case of 
SN 1006. Finally, radio, X-ray and y-ray images are synthesized 
at 1 GHz, 3 keV, and 1 TeV, respectively. It is worth to empha- 
size that all the above parameters are not expected to influence 
the pattern of the asymmetries induced by a non-uniform ISMF 
on which the present paper is focused (see, also, Sect. 14.61 for a 
discussion on the influence of £max,|| on the remnant morphol- 
ogy). 

4.2. Asymmetries in the remnant morphology: the reference 
case 

In Paper I, we analyzed the asymmetries induced by a non- 
uniform ISMF in the radio morphology of the remnant. In par- 



6 Note that Ef t \\ depends on the magnetic field strength at paral- 
lel shock which is different in the three configurations of unperturbed 
ISMF explored here due to the magnetic field gradient (in all the cases, 
we fix the magnetic field strength at the center of the SN explosion). 
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Figure 3. Azimuthal profiles of the maximum energy £ max com- 
puted in run Grad-BZ-gl (solid line) and Grad-BX-gl (dashed 
line) when the aspect angle is 0b = 90° and the gradient of mag- 
netic field strength lies in the plane of the sky (0vb = 0° in run 
Grad-BZ-gl and 0vb = 90° in run Grad-BX-gl). The shock is 
parallel around 90° and 270° and perpendicular around 1 80° and 
360°. In both models, the contrast C ma x > 1. 



ticular, we found there that asymmetric BSNRs are produced 
if a gradient of the ambient magnetic field strength IB is not 
aligned with the LoS. In this section we extend our analysis to 
non-thermal X-rays and IC y-rays. To this end, we synthesize the 
synchrotron and IC emission, considering each of the three cases 
of variation of electron injection efficiency with shock obliquity 
(quasi-perpendicular, isotropic, and quasi-parallel particle injec- 
tion). Also we assume the adiabatic index to be y = 5/3; the ef- 
fects of lower y values on the remnant morphology are explored 
in Sect. 14.31 Here the maximum energy of electrons is calcu- 
lated at each point as £ max = min[£ max ,i, £ max ,2, E max3 ] (where 
indexes 1,2,3 correspond respectively to loss-limited, time- 
limited and escape-limited models; see discussion in Sect. 13. lb . 
For the set of parameters chosen for our simulations, it turns 
out that the loss-limited model is dominant at all obliquity an- 
gles, thus simplifying the analysis of non-thermal images in this 
section. Figure [3] shows £" max versus the azimuthal angle (the 
azimuth is measured counterclockwise from the "north" of the 
remnant) for runs Grad-BZ-gl and Grad-BX-gl. In both cases, 
E max is characterized by two maxima where the ISMF is paral- 
lel to the shock normal (around 90° and 270°); thus the contrast 
of E mSLX is Cmax > 1 • The strength of the unperturbed ISMF is 
the largest at 180° (90°) and the lowest at 360° (270°) in run 
Grad-BZ-gl (Grad-BX-gl) due to the magnetic field gradient. 
The latter determines the asymmetries in the azimuthal profile 
of E max : in run Grad-BZ-gl, the two maxima are converging on 
the side where the field is the most intense, the gradient being 
perpendicular to the average magnetic field; in run Grad-BX-gl, 
the two maxima have different intensities (with the largest where 
the magnetic field is the lowesQ), the gradient being parallel to 
the average magnetic field. 

As an example, Figs. [Hand [5] show the maps of synchrotron 
radio, X-ray, and IC y-ray surface brightness at t = 1 kyr, in each 
of the three injection models (quasi-perpendicular, isotropic, and 
quasi-parallel). The aspect angle is 0b = 90° in all images, i.e. 
the ambient magnetic field is perpendicular to the LoS; the angle 
0vb is 0° for run Grad-BZ-g 1 and 90° for Grad-BX-g 1 . 



In the loss-limited model E max oc B 1/2 , see Sect. 13.11 
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Figure 4. Maps of synchrotron radio (left), X-ray (center), and 
IC y-ray (right) surface brightness (normalized to the maximum 
of each map) at t = 1 kyr synthesized from run Grad-BZ-gl, 
assuming randomized internal magnetic field. The relevant an- 
gles are 0b = 90° and 0vb = 0°. The figure shows the quasi- 
perpendicular (top), isotropic (middle), and quasi-parallel (bot- 
tom) particle injection models. The adiabatic index is y = 5/3. 
The average ambient magnetic field is along the horizontal axis; 
the gradient of magnetic field strength is along the vertical axis. 

The main factors affecting the azimuthal variations of surface 
brightness are the variations of: injection efficiency f g (S ) an d 
magnetic field B s (® ) in the radio band; f g (& ) 9 B s (® ) and max- 
imum energy E max (® ) in the X-ray band; f g (S ) and E max (® ) 
in the IC y-ray band. Therefore, the morphology of the rem- 
nant in the three bands can differ considerably in appearance. In 
the radio and in the X-ray band, the remnant shows two lobes 
located at perpendicular shocks in the quasi-perpendicular and 
isotropic models (i.e. where the magnetic field is larger), and 
at parallel shocks in the quasi-parallel model (i.e. where emit- 
ting electrons reside). The lobes are much thinner in X-rays 
than in radio because of the large radiative losses at the high- 
est energies that make the X-ray emission dominated by radii 
closest to the shock. In the y-ray band, the remnant morphol- 
ogy changes significantly in the three injection models: it is al- 
most ring-like (with two faint minima at parallel shocks) when 
the injection is quasi-perpendicular; the morphology shows two 
lobes located at parallel shocks when the injection is isotropic, 
at variance with the lobes in radio and X-rays that are located 
at perpendicular- shocks (i.e. bright y-ray lobes correspond to 
dark radio and X-ray areas); the morphology is characterized 
by two narrow bright lobes almost superimposed to those in 
radio and X-rays when the injection is quasi-parallel. A ring- 
like y-ray morphology is compati ble with those found b y HESS 
in the SNRs RX J171 3.7-3946 (lAharonian eTai1 l2006h and RX 
J0852.0-4622 (Vela Jr.: lAharonian et al.l2007bl) where y-rays are 
detected virtually throughout the whole remnant and the emis- 
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Figure 5. As in Fig.|4]for run Grad-BX-gl . Both the average am- 
bient magnetic field and the gradient of magnetic field strength 
are along the horizontal axis. The relevant angles are cpB = 90° 
and 0vb = 90°. 



sion is found to resemble a shell structure. On the other hand, 
the bipolar y-ray m orphology of SN 1006 revealed by HESS 
(lAcero et al.l [2010), with the bright lobes strongly correlated 
with non-thermal X-rays, may be easily reproduced in the polar- 
caps scenario (quasi-parallel injection). 

The effects of the non-uniform ISMF on the remnant mor- 
phology in the X-ray band are similar to those discussed in 
Paper I for the radio band: remnants with two non-thermal X-ray 
lobes of different brightness (upper left panel in Fig.|4]and lower 
left panel in Fig.0 are produced if a gradient of ambient mag- 
netic field strength is perpendicular to the lobes; remnants with 
converging similar non-thermal X-ray lobes (lower left panel in 
Fig. |4] and upper right panel in Fig. 0) are produced if the gra- 
dient runs between the two lobes. Analogous asymmetries are 
found in the y-ray morphology of the remnant although the de- 
gree of asymmetry is less evident. Note however that, in the case 
of isotropic injection, the y-ray lobes are converging on one side 
when radio and X-ray lobes are characterized by different bright- 
ness (see Fig.©. Th is is the consequenc e of the "limb-inverse" 
property in y-ray dPetruk et al. 2009a). In general, this prop- 
erty is valid not only in the case of isotropic injection; this type 
of injection is just the more prominent case. In fact, the critical 
quantities determining the "limb-inverse" property are the con- 
trasts between electron injection, ISMF, and model of E max . For 



8 The "limb-inverse" property in y-rays is determined for isotropic 
injection because the magnetic field affects the downstream distribu- 
tion of IC y-ray emitting electrons which i s steeper where the magnetic 
field is stronger. The reader is referred to lPetruk et al.1 (l2009ah for more 
details. 
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Figure 6. Azimuthal profiles of the synchrotron 
radio (red), X-ray (blue), and IC y-ray (green) 
surface brightness synthesized from runs Grad- 
BZ-gl (left; the relevant angles are B = 
90° and VB = 0°) and Grad-BX-gl (right; 
B = 90° and 0vb = 90°), assuming quasi- 
perpendicular (top), isotropic (middle), and 
quasi-parallel (bottom) injection models. The 
adiabatic index is y - 5/3. The azimuth is 
measured counterclockwise from the north (see 
Figs. |4] and [5]). The corresponding azimuthal 
profiles of E max used to derive the curves of this 
figure are in Fig. [3] 



instance, in the case of uniform ISMF, the azimuthal contrast in 
IC y-ray brightness is roughly 



S || injection] | 

~~~ OC - ; ; exp 

S± injection^ 



1 



1 



(21) 



injection] | 
injection. 



■ exp 



where E m is the electron energy which gives the maximum con- 
tribution to IC emission at a considered frequency and sub- 
scripts refer to positions along the limb where the ambient 
magnetic field is either parallel (||) or perpendicular (_L) to 
the shock normal. Even in the case of quasi-parallel injection 
(injection||/injection ± > 1), the contrast S\\/S± depends on the 
contrast of E mSiX \ the ratio C ma x = £max,||/£max,± may lead to an 
exponential term either > 1 or < 1, leading to S\\/S± either > 1 
or < 1. 

Another interesting feature characterizing the IC y-ray mor- 
phology of the remnant is the inversion of the asymmetry when 
the two lobes have different brightness (i.e. a gradient of mag- 
netic field strength is perpendicular to the lobes). This feature is 
evident in the upper panels of Fig. |4] and in the lower panels of 
Fig. [5j the brightest y-ray lobe is located where both the radio 
and the X-ray lobes are fainter. As discussed in detail below in 
Sect. 14.41 this is due to the fact that, in the synthetic images pre- 
sented in this section, E max depends inversely on the pre-shock 
ambient magnetic field strength (see Eq. [TO] and Fig. [3]) and its 
contrast is C ma x > 1 • 

Figure [6] shows the azimuthal profiles of the synchrotron ra- 
dio, X-ray, and IC y-ray surface brightness synthesized from 



runs Grad-BZ-gl and Grad-BX-gl for the three injection mod- 
els when the relevant angles are 0b = 90° and 0vb = 0° 
for run Grad-BZ-gl and V b = 90° for Grad-BX-gl. In the 
quasi-parallel scenario, the non-thermal lobes are rather nar- 
row azimuthally. Note the "limb-inv erse" property in y-r ay s for 
isotropic injection as discussed by Petr uk et aD (l2009al) . Note 
also the "asymmetry-inverse" property in y-rays when the two 
lobes have different brightness. In general we find that the de- 
gree of asymmetry (whatever the pattern of asymmetry - either 
different brightness or convergence of the lobes - is) induced by 
VB in the remnant morphology is different in the three bands 
(see Sect. 14.51 for a discussion on the degree of asymmetry of 
the remnant in the different bands); in particular, the IC y-ray 
emission appears to be the less sensitive to the gradient. 



Useful parameters to quantify the degree of asymmetry of 
the remnant are those defined in Paper I: the azimuthal intensity 
ratio R max > 1 , i.e. the ratio of the maxima of intensity of the two 
lobes as derived from the azimuthal intensity profiles (a measure 
of different brightness of the lobes; R mSiX > 1 in case of asym- 
metry), and the azimuthal distance #d, i.e. the distance in deg 
of the two maxima (a measure of the convergence of the lobes; 
#d < 180° in case of asymmetry). For instance, in the case of 
quasi-parallel injection in Fig. [6] (lower panels), we find that the 
azimuthal distance #d ranges from 148° in y-rays and X-rays to 
134° in radio for run Grad-BZ-gl, and the azimuthal intensity 
ratio R max ranges from 1.4 in y-rays to 1.8 in radio and X-rays 
for run Grad-BX-gl. 
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Figure 7. As in Fig. @] for runs Grad-BZ-gl (7 = 5/3), Grad-BZ-g2 (7 = 4/3), and Grad-BZ-g3 (7 = 1.1). The figure shows 
synchrotron radio (left panels), X-ray (center panels) and IC 7-ray (right panels), assuming either quasi-perpendicular (upper panels) 
or quasi-parallel (lower panels) injection models. Each panel shows only one half of the remnant which is symmetric with respect 
to the vertical axis. 



4.3. Dependence on the adiabatic index 

Petruk et al. (2010b, submitted to MNRAS) analyzed the effect 
of 7 on non-thermal images of SNR expanding through homo- 
geneous ISM and uniform ISMF. They showed that reducing the 
value of 7, the synchrotron brightness of the remnant is modi- 
fied by increased radiative losses of emitting electrons, due to 
increased compression of B, which results in thinner radial pro- 
files of brightness. Figure [7] shows maps of synchrotron radio, 
X-ray, and IC 7-ray emission for the case of ISMF character- 
ized by a gradient of field strength perpendicular to the average 
magnetic field and different values of the adiabatic index 7 (runs 
Grad-BZ-gl, Grad-BZ-g2, and Grad-BZ-g3). As expected, the 
index 7 determines both the shock compression ratio cr and the 
distance of the contact discontinuity from the blast wave position 
Z) C( j (see also Fig [T] in the case of uniform ISMF): the smaller 7, 
the larger cr (and the larger the radiative losses of emitting elec- 
trons) and the shorter D c $. As shown in the figure, the main effect 
of smaller 7 is to make thinner the lobes emitting synchrotron 
emission in the three bands. In particular, in the extreme case of 
7 = 1 . 1 , the lobes are so thin that they are largely perturbed by 
the hydrodynamic instabilities forming at the contact disconti- 
nuity, the typical size of the instabilities being comparable with 
Z) C( j. The adiabatic index slightly influences also the azimuthal 
thickness of the lobes, especially in the quasi-parallel case: the 
smaller 7, the narrower this thickness. Nevertheless, the adia- 
batic index does not change significantly neither the degree nor 
the pattern of asymmetry of the remnant morphology caused by 
the gradient of magnetic field strength. 

4.4. Dependence on the maximum energy 

In Sect. 14.21 we have presented the remnant morphology for 
a reference case for which E max calculated at each point as 



£max = mm[£ max j,£max,2,£max,3] (where indexes 1,2,3 cor- 
respond respectively to loss-limited, time-limited and escape- 
limited models; see discussion in Sect. 13. II) is characterized by a 
contrast C m ax > 1 • Here we generalize our study by considering 
some arbitrary smooth variations of £ max versus obliquity with 
the goal to see how different trends in the obliquity dependence 
of ^max influence the visible morphology of SNRs. That is, we 
do not use here the particular prescriptions for £ max from Sect. 
13.11 but simply assume that the acceleration physics is able to 
operate to produce an £ max of prescribed properties. In fact, the 
model of £" m a X and its obliquity dependence may affect both the 
degree and the pattern of asymmetry of the remnant morphol- 
ogy. A critical point is if £ max depends directly or inversely on 
the magnetic field strength. Moreover the remnant morphology 
in the various bands can be characterized by different features 
if the contrast C ma x is either > 1 or < 1. As an example, Fig. [8] 
shows the azimuthal profiles of three arbitrary models of E max 
characterized by different dependencies on the obliquity angle: 
in model A, E m2LX °c B and its contrast is C ma x < 1 (solid line); 
in model B, E max oc B~ 112 and C ma x > 1 (dotted); in model C, 
£max 00 B and C ma x is < 1 on the side of the remnant with the 
strongest ISMF strength and > 1 on the side with the lowest 
field strength. Note that model B coincides with the model of 
£max computed for the reference case in Sect. 14.21 (see Fig. [3]). 
The asymmetries on the profiles of E mSLX are introduced by the 
gradient of ambient magnetic field, as discussed in Sect. 14.21 

From the models of E max shown in Fig. [8l we synthesized 
maps of synchrotron radio, X-ray, and IC 7-ray emission. We 
find that the remnant morphology is very sensitive to the model 
of E mdLX when a gradient of magnetic field strength is perpen- 
dicular to the lobes and the latter are characterized by different 
brightness. The asymmetry between the two lobes can be re- 
duced in the X-ray band or even inverted in the IC 7-ray band 
when £ max depends inversely on the pre- shock ambient mag- 
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Figure 8. As in Fig. [3] but for three different models of E max 
characterized by smooth variations with the obliquity angle. The 
azimuth is measured counterclockwise from the "north" of the 
remnant. The shock is parallel around 90° and 270° and per- 
pendicular around 180° and 360°. In model A (solid line), the 
contrast of E max is C m ax < 1 , in model B (dotted) C m ax > 1 , and 
in model C (dashed) C ma x is < 1 on the side of the remnant with 
the strongest ISMF strength and > 1 on the side with the lowest 
field strength (see text). 



netic field strength, namely in the case of model B in Fig. [8] In 
particular, this model of E max leads to the "asymmetry-inverse" 
property in y-rays already discussed in Sect. 14.21 for our refer- 
ence case (see upper right panel in Fig. |4] and lower right panel 
in0. 

Fig. [9] shows the azimuthal profiles of the IC y-ray surface 
brightness synthesized from runs Grad-BZ-gl and Grad-BX-gl 
when the lobes have different brightness, for the three models 
of E max reported in Fig. [8] In the case of the IC surface bright- 
ness, the asymmetry-inverse property is evident when E mSLX is 
the largest where the magnetic field strength is the lowest (see 
model B in Fig. [8]). This is due to the fact that the IC emissiv- 
ity i(e) weakly depends on B. In the case of the non-thermal 
X-ray surface brightness, the inverse dependence of E max on B 
partially contrast the dependence of the non-thermal X-ray i(e) 
on B, reducing the degree of asymmetry between the lobes. It is 
worth to note that, if E max is high enough in regions with weak 
magnetic field, than the inversion of asymmetry may be present 
even in the X-ray band. 
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Figure 9. Azimuthal profiles of the synchrotron radio (red lines) 
and IC y-ray surface brightness (green) synthesized from run 
Grad-BZ-gl with quasi-perpendicular injection (top), and from 
run Grad-BX-gl with quasi-parallel injection (bottom); the as- 
pect angle is 0b = 90°. The models of E max adopted to synthe- 
size the non-thermal emission are those shown in Fig. [8j model 
A (solid line), B (dotted) and C (dashed). 



On the other hand, we also found that when the non-uniform 
ISMF leads to non-thermal lobes converging on one side (i.e. 
when a gradient of ISMF is running between the lobes) the 
model of E max does not affect significantly the degree and the 
pattern of asymmetry of the remnant morphology (see lower 
panel in Fig. [9]). 



4.5. Dependence on the orientation of ISMF gradient 

As expected, the degree of asymmetry of the remnant morphol- 
ogy depends on the orientation of V# with respect to the plane 
of the sky. In the case of run Grad-BZ-gl, Fig. [TOl shows the az- 
imuthal intensity ratio R mSLX and the azimuthal distance #d vs. 
the angle 0vb, for an aspect angle 0b = 90° and for differ- 
ent trends in the obliquity dependence of E mSLX (exploring the 
contrasts C ma x either > 1 or < 1; see Fig. [8]). The asymme- 
tries are the largest when VB lies in the plane of the sky (i.e. 
0vb = 0°), whereas no asymmetries are present when VB is 
along the LoS (i.e. 0vb = 90°). In all the intermediate cases, 
the degree of asymmetry is determined by the component of VB 
lying in the plane of the sky. Note that the remnant morphology 
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Figure 10. Azimuthal intensity ratio R m2LX (i.e. the ratio of the 
maxima of intensity of the two lobes around the shell) and az- 
imuthal distance #d (i.e. the distance in deg of the two maxima 
of intensity around the shell) vs. the angle between VB and the 
vertical line passing through the remnant center 0vb, for an as- 
pect angle 0b = 90°, and for the three models of E max shown 
in Fig. [HJ The run is Grad-BZ-gl. For isotropic injection, curves 
for R max in y-rays and curves for #d in radio and X-rays are not 
shown, the values being R max - 1 and #d = 180°, respectively at 
all 0vB- 



shows only one kind of asymmetry when the injection is quasi- 
perpendicular or quasi-parallel and the aspect angle is 0b = 90°. 
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Figure 11. Maps of synchrotron radio (left), X-ray (center), and 
IC y-ray (right) surface brightness synthesized from run Grad- 
BZ-gl, assuming quasi-perpendicular (top), and quasi-parallel 
(bottom) injection models. The adiabatic index is y = 4/3. The 
maps have been synthesized adopting model B of E mSLX shown 
in Fig. [8j The relevant angles are 0b = 45° and 0vb = 45°. The 
angle between (B) and V|(B)| is 90°. 



On the other hand, the lobes have different brightness in radio 
and non-thermal X-rays and are converging in IC y-rays when 
the injection is isotropic due to the "limb-inverse" property. 

In general we find that the degree of asymmetry (whatever 
the pattern of asymmetry - either different brightness or conver- 
gence of the lobes - is) induced by VB in the remnant morphol- 
ogy is different in the three bands: the non-thermal X-ray (IC 
y-ray) emission appears to be the most (less) sensitive to the gra- 
dient. This happens because the emissivity /(e) depends directly 
on the magnetic field strength (see Eq.0 only in the synchrotron 
emission process (no in the IC process). Consequently, the IC y- 
ray emission shows a weaker dependence on the VB. In fact, IC 
brightness depends on B indirectly, through radiative losses of 
electrons: larger B induces decrease of the number of electrons 
emitting IC y-ray s. Note that the sensitivity on VB depends also 
on the energy of photons, and on the reduced fiducial energy 
Ef, which is the measure of efficiency of the role of radiative 
losses in modification of the downstream evolution of emitting 
electrons. Note also that the degree of asymmetry of the remnant 
morphology can be significantly reduced when the E max contrast 
is Cmax > 1 (e.g. model B in Fig. [8]). In fact, the asymmetry 
reduction is mainly due to the dependency of E mSLX on the mag- 
netic field strength. As discussed in Sect. 14.41 the asymmetries in 
the remnant morphology can be reduced or even inverted when 
E max depends inversely on the pre- shock ambient magnetic field 
strength (which is the case in model B). 

When the VB is not aligned with the average ambient mag- 
netic field (for instance in the case of run Grad-BZ-gl), the pro- 
jection of the VB in the plane of the sky has (for generic val- 
ues of 0b and 0vb) a component perpendicular to the projected 
lobes and one running between them. In this case both kind of 
asymmetries (lobes converging on one side and with different 
brightness) are expected in the remnant morphology. As an ex- 
ample, Fig. [TT] shows the synchrotron radio, X-ray, and IC y-ray 
images synthesized from run Grad-BZ-gl, for different injection 
models. The relevant angles are 0b = 45° and 0vb = 45°. 
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Figure 12. As in Fig. for the azimuthal profiles of the syn- 
chrotron X-ray (blue) and IC y-ray (green) surface brightness 
synthesized from run Grad-BZ-gl, for £ max ,|| = 26 TeV (solid), 
5 TeV (dotted), and 1 TeV (dashed). The surface brightness is 
synthesized deriving the maximum energy £ max as described in 
Sect. 14. 2t its azimuthal profile is similar to that shown in Fig. [3] 
(solid line) but with a different value of £ m ax,n- In all the cases, 
£f,n = 12 TeV. 



4.6. Dependence on the value of E, 



max, 1 1 



The calculations presented above assume the free parameter 
£max,n = 26 TeV larger than the fiducial energy Ef t \\ = 12 TeV, 



implying that the electron energy losses are mainly due to ra- 
diative losses (see Sect. 13.2b . The variation of the energy distri- 
bution N(E, a, t) of electrons in Eq. [19] is influenced by radia- 
tive losses of electrons that are, therefore, important for the sur- 
face brightness distribution of the remnant in X-rays and y-rays. 
Consequently, the choice of E max ,\\ ma Y influence both the de- 
gree and the pattern of asymmetry of the remnant morphology. 
In particular, in cases with £ max < £f, we expect that £ ra d — > 1 
(i.e. electrons have negligible radiative losses) and the electron 
energy losses are mainly due to adiabatic expansion. This issue 
is investigated by considering the reference case discussed in 
Sect. 14. 21 and two additional cases for which £ max < £>, namely 
£max,n = 5 TeV and is max, y = 1 TeV. Figure [T21 shows the az- 
imuthal profiles of the synchrotron X-ray and IC y-ray surface 
brightness synthesized from run Grad-BZ-gl, for these values of 
£max,n together with the case with £ max ,n = 26 TeV (the reference 
case). Note that E max is calculated at each point of the domain as 
described in Sect. l4.2l but for different values of £max,|| (see also 
Sect. 13.11) . The figure shows that for decreasing values of E max ,\\, 
the contrast of emission increases, the effect being the largest for 
IC y-ray emission than for synchrotron X-rays. Nevertheless, the 
degree and the pattern of asymmetry of the remnant morphology 
induced by the gradient of ISMF are only slightly influenced by 
the value of E m *x ii. 



5. Summary and conclusions 

We developed a numerical code (remlight) to synthesize the 
synchrotron radio, X-ray, and IC y-ray emission from MHD sim- 
ulations, in the general case of a remnant expanding through a 
non-uniform ISM and/or a non-uniform ISMF. As a first appli- 
cation of remlight, we coupled the synthesis code to the MHD 
model discussed in Paper I (extended to include an approximate 
treatment of upstream magnetic field amplification and the ef- 
fect of shock modification due to back reaction of accelerated 
CRs) and investigated the effects of a non-uniform ISMF on the 
remnant morphology in the X-ray and y-ray bands. Our findings 
lead to several conclusions: 

- A gradient of ISMF strength induces asymmetries in both the 
X-ray and y-ray morphology of the remnant if the gradient 
has a component perpendicular to the LoS. In general, the 
asymmetries are analogous to those found in Paper I in the 
radio band, independently from the models of electron in- 
jection and of maximum energy of electrons accelerated by 
the shock. In the y-ray band, the asymmetry in the remnant 
morphology is inverted with respect to those in the radio and 
X-ray bands if the model of £ max depends inversely on the 
pre- shock magnetic field strength and its contrast is C ma x > 1 
(e.g. model B in Fig. [8]): the brightest y-ray lobe is located 
where both the radio and the X-ray lobes are the faintest. 

- The non-thermal lobes are characterized by different bright- 
ness when a gradient of ISMF strength is perpendicular to 
the lobes; they are converging on one side when a gradient 
of ISMF is running between them. In the general case of a 
gradient with components parallel and perpendicular to the 
lobes, both kinds of asymmetry may characterize the rem- 
nant morphology. 

- The non-thermal X-ray emission is confined in very thin 
limbs because of the large radiative losses at high energy 
and, in general, is the most sensitive to non-uniform ISMF. 
In fact the remnant morphology in this band shows the high- 
est degree of asymmetry among the images synthesized in 
the three bands of interest (i.e. radio, X-ray, and y-ray), ex- 
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cept when E max depends inversely on the pre- shock magnetic 
field strength. In the latter case, the asymmetries in the X-ray 
band can be significantly reduced. 

The IC y-ray emission is weakly sensitive to the non-uniform 
ISMF, the degree of asymmetry being the lowest in the 
three bands considered. The remnant morphology is al- 
most ring-like for quasi-perpendicular injection, shows the 
"limb-inverse" property discussed by iPetruk et aD (l2009ah 
for isotropic injection (i.e. bright y-ray lobes correspond 
to dark radio and X-ray areas), and is bilateral for quasi- 
parallel injection. The "limb-inverse" property implies, for 
instance, that y-ray lobes are symmetric and converging on 
one side when radio and X-ray lobes have different bright- 
ness (see Fig. [4]). In case E max depends inversely on the 
pre-shock ambient magnetic field strength, the asymmetries 
in the IC y-ray morphology can be inverted; for instance, 
brightest y-ray lobes can be located where both radio and 
X-ray lobes are fainter. Note that the y-ray morpho logy 
of the SNRs RX J 1713.7-3946 (lAharonian et alJ l2006) and 
RX J0852.0-4622 (lAharonian et alJl2007bb could be repro- 
duced in the equatorial-belt scenario (the injection is either 
quasi-perpendicul ar or isotrop ic), whereas the morphology 
of SN 1006 (lAcero et al.ll2010h is compatible with that pre- 
dicted in the polar-caps scenario (quasi-parallel injection). 



value of an "effective" y is provided in each point of the spatial 
domain^. 

Note also that the MHD model adopted here follows the evo- 
lution of the remnant during the adiabatic phase and, therefore, 
its applicability is limited to this evolutionary stage. In the radia- 
tive phase, the high degree of compression suggested by radia- 
tive shocks leads to increase in the synchrotron emission bright- 
ness due to compression of ambient magnetic field and electrons. 
Since our model neglects the radiative cooling of the shocked 
gas, it is limited to compression ratios derived from y and, there- 
fore, it is not able to simulate this mechanism of limb bright- 
ening. Nevertheless, the model is appropriate to describe young 
SNRs that are those from which non-thermal emission is com- 
monly detected. 
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Note that, although the MHD model presented here does not 
include self-consistently shock modification and magnetic field 
amplification, we adopted an approximate treatment of both pro- 
cesses. Magnetic field amplification could result from stream- 
ing instability excited by the accelerated particles upstream of 
the shock or, alternatively, the magnetic fields could be ampli- 
fied in a purely hydrodyna mic way in the downstream plasma 
(iGiacalone & Jokipiill2007l) . In both cases, the shock is expected 
to be modified due to the dynam ical reaction of the amplified 
magnetic field (see, for instance. [Ferrand et al.ll2010l for an hy- 
drodynamic model including back-reaction of accelerated CRs). 
In this paper, we approach the effect of shock modification by 
considering different values of the adiabatic index y (namely, 
5/3, 4/3, 1.1) and the effect of upstream magnetic field ampli- 
fication by considering the ambient magnetic field strength en- 
hanced by x 10 in the neighborhoods of the remnant (the unper- 
turbed field strength commonly expected is a few fiG). The main 
effect of y is to change the compression ratio of the shock and the 
distance of the contact discontinuity from the blast wave posi- 
tion. In the simplest case considered here, namely the modifica- 
tion on y and the upstream magnetic field amplification are both 
isotropic with no dependence on the obliquity angle, we found 
that the modified y and the amplified field influence mainly the 
absolute values of non-thermal emission but not the large scale 
morphology of the remnant and the pattern of asymmetries in- 
duced by a non-uniform ISMF. The results presented here, there- 
fore, are only valid in this case. Conversely, we expect a signif- 
icant effect of the modified y as well as of the amplified field 
on the remnant morphology if the shock modification and/or 
upstream magnetic field amplification depend on the obliquity. 
This issue deserves further investigation in future studies. 

It is worth to emphasize that the calculations provided in this 
paper (and implemented in the remlight code) to synthesize the 
non-thermal emission from MHD simulations consider a generic 
adiabatic index y. The remlight code therefore can be easily 
coupled with a model including the back-reaction of accelerated 
CRs and synthesize the non-thermal emission consistently if the 
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Now, the relation r(a,R ) is approximatecQ using the method described by 
Hnatvk & Petruk (1999): 



r(a,R) (a\r 2 3 4 

=1 — 1 (1 + a\v + a^v + a^v + a^v 

R \R/ 



(A.7) 



where v = (R - a)/R and if/ = (y - l)/y. The parameters a\, «3, and are 
expressed as: 



a\ = -r a , s + if/ , 



a 2 = - (Rr aa , s - 2if/r a , s + if/(i// + 1)) 



Cl2 



\ (~R 2 r aaa ,, + 3i/tf r aa , s - 3<A0A + l)r a , s + + 1)(<A + 2)) , 



«4 = C - (1 + a\ + <22 + a^) , 



(A.8) 
(A.9) 

(A. 10) 
(A. 11) 



where C reflects the variation of r(a) around the center of the SNR. We adopt 
C = Ca where Ca is given by the self-similar S edov solution for a spherical 
shock (for details see Sect. 4.3 and Appendix in Petruk 2000 and references 
therein): 



Appendix A: Calculation of integral in Eq.FRl 

The function I(a, t) in Eq.[l4]is expressed as ( Revnoldsl l 19981) 



I(a, t) = (j\ 



r f B eff (a,t') 2 ( p(a,f) \ 1/3 df 
X #eff,s(0 2 \p(a,t)J t 



(A.l) 



where B 2 S = B 2 + #£ MB is the "effective" magnetic field introduced to account 
for the energy losses of electrons due to IC scatterings on the photons of CMB. 



The integral JA.U is rather CPU consuming because it requires to know, 
with high enough time resolution, the history of each parcel of gas inside the 
SNR since its shocking time. To reduce the computational cost, we calculate it 
approximately, changing integration on df to dR' = V^(f)df ', where R and V s h 
are the shock position and velocity, respectively, and using some MHD properties 
of the fluid. 

We calculate I(a, t) using an analytic description of mass density and mag- 
netic field evolution inside the SNR which expands through a non-uniform ISM 
and/or ISMF. The continuity equation p (a)a 2 da = p(a)r 2 dr results in 



p(a, t) = p Q (a) [^—^ j r a(«, t) 



(A.2) 



wher e r a (a, t) is the derivative of r(a, t) with respect to a; the density term in Eq. 
lAJlis 



C A - 



y + 1 



P(0) 



-i/r 



1/3 



(A. 12) 



P(0) is the plasma pressure at the center of the remnant divided by its post-shock 
value 



^(0) 



= / 1 \ 6/5 /r+i\ 6/5 " y/(2 " y) / (2 r+ i)(y+D \ ( - 

\2) \ 7 J \ 7(7-7) J 



■2+5/(2-r)K 



( = 



7+1 



3( r -l) + 2 



2 7-1 
' 5 + 2 r +l 



(A. 13) 



(A. 14) 



Thus, we derive C A (7 = 5/3) = 1.083, C A (7 = 4/3) = 1.055 and C A (y = 1.1) = 
1.021. The expressions for the derivatives r a;S , r aajS , r aaa;S in Eqs. |A.8|A.10| as 
functio ns of R, R, R and R^ are given in Appendix A2 of Hnatvk & Petruk 
(H999h . 

Finally, w e calculate V s h(^) in Eq. lA.ll as well as R and R^ 3 \ using the 
Hnatvk (1987) approximate analytical formula for the strong shock in a non- 
uniform medium (see also Sect. 2.1 in Hnatvk & Petruk 1999). 

The integral I(a, t) can be calculated rather simply in the case of a SNR 
expanding through uniform ISM and ISMF. We therefore test our calculation 
of I(a,t) by comparing the approximate values derived from Eq. I A. 6 1 wit h the 
exact ones derived from the Sedov solution in the case of 7 = 5/3. Figure IaH 
compares the exact and approximate values of I (a, t) in the limits of parallel and 
perpendicular shocks. Note that the approximate values are very accurate at radii 
close to the shock front where most of the non-thermal emission originates. 



p(a,t') = r(a,t) 2 r a (a,t) 
p(a,t) r(a,f) 2 r a (a,f) 



(A3) 



The magnetic field in Eq. lA.ll can be expressed as B(a, t) 2 = B\\(a, t) 2 +B ± (a, t) 2 , 
where B\\ and B ± are the components of magnetic field parallel and perpendicular 
to the shock normal, respectively. These two components follow the magnetic 
flux conservation B\\dos = const, where dos is a surface element, and the flux- 
frozen condition B ± (r)rdr = const: 



B\\(a,t) = B lo (a) 



r(a, t) 2 



(A.4) 



B±(a , t) = B±Aa) ^e^ = B±Aa) « 



a p (a) 



r(a, t)r a (a, t) 



(A.5) 



Thus, the magnetic field and the mass density in Eq. IA.1I can be expressed 
through the relation r(a, t) between Eulerian and Lagrangian coordinates of a 
parcel of gas and its derivative, r a (a, t). Considering that r(a, t) and r a (a, t) can 
be expressed in terms of t he dy namical characteristics of the shock (i.e. as r(a, R) 
and r a (a,R)), the integral [Pl may be calculated as follows: 



I(a,t). 



o\ r R B(a,R') 2 I r(a,R) 2 r a (a,R) \ 1/3 
t J R . B s (R) 2 \r(a,R') 2 r a (a,R>)) 



dR' 



(A.6) 



10 The approximation dA.7l ) is developed to give exact values of 
derivatives up to the third order at the shock and to the first order at 
the center. 



16 



S. Orlando et al.: Effects of non-uniform ISMF on X-ray and y-ray morphology of SNRs 




0.0 0.2 0.4 0.6 0.8 
x/r SNR (parallel shock) 



0.0 0.2 0.4 0.6 0.8 
z/r SNR (perpendicular shock) 



1.0 




1.0 



Figure A.l. Self-similar approximate and exact radial profiles 
of the integral 1(a) when the ambient magnetic field is either 
parallel (upper panel) or perpendicular (lower panel) to the shock 
normal, and y = 5/3. 



